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Abstract 



We carry out a systematic investigation of twist-two spin dependent structure functions 
of the nucleon within the framework of the chiral quark soliton model (CQSM) by paying 
special attention to the role of chiral symmetry of QCD. The importance of chiral symmetry 
is illustrated through the good reproduction of the recent SLAC data for the neutron spin 
structure function gi{x, Q^). We also observe substantial difference between the predictions of 
the longitudinally polarized distribution functions and those of the transversity distribution 
functions. That the chiral symmetry may be responsible for this difference is seen in the 
isospin dependence of the corresponding first moments, i.e. the axial and tensor charges. 
The CQSM predicts / — 0.25 for the ratio of the isoscalar to isovector axial charges, 
while g^^ /g^^ — 0.46 for the ratio of the isoscalar to isovector tensor charges, which should 
be compared with the prediction g'^^ IdA ~ 9t'^ /dr^ = 3/5 of the constituent quark model or 
of the naive MIT bag model without proper account of chiral symmetry. Another prominent 
prediction of the CQSM is the opposite polarization of the u and d antiquarks, thereby 
indicating the SU(2) asymmetric sea quark (spin) polarization in the nucleon. 
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1 Introduction 



Undoubtedly, the so-called "nucleon spin crisis" caused by the EMC measurement in 1988 
is one of the most exciting topics in the field of hadron physics [1] . The recent renaissance of 
nucleon structure function physics is greatly owing to this epoch-making finding. Naturally, 
the physics of nucleon structure functions has two different aspects. One is a perturbative 
aspect, while the other is a nonperturbative aspect. Because of the asymptotic freedom of 
QCD, the Q^-evolution of quark distribution functions can be controlled by the perturbative 
QCD at least for large cnoTigh Q"^ [2]. However, the perturbative QCD is entirely powerless 
for predicting distribution functions themselves. Here we need to solve nonperturbative QCD 
in some way. Unfortunately, we have no reliable analytical method for handling this aspect 
of QCD. For the present moment, we are then left with two tentative choices. One is to rely 
upon lattice QCD, while the other is to use effective models of QCD. If one takes the first 
choice, one must first evaluate infinite towers of moments of distribution functions, since the 
direct calculation of distribution functions does not match this numerical simulation method 
[3]. Here we take the second choice, which allows us a direct calculation of quark distribution 
functions. Still, there are quite a lot of effective model of baryons. We advocate that the 
chiral quark soliton model (CQSM) is a unique model of baryons which has several appealing 
features not possessed by other models of baryons, especially when applied to the physics 
of quark distribution functions. First of all, it is an effective model of baryons maximally 
incorporating spontaneous chiral symmetry breaking of QCD vacuum [4-6]. The nucleon in 
this model is a composite of three valence quarks and infinitely many Dirac sea quarks moving 
in a slowly rotating M.F. of hedgehog shape. As a natural consequence, it automatically 
simulates cloud of pions surrounding the core of three valence quarks. Nevertheless, since 
everything is described in terms of effective quark fields only, we need not worry about a 
double counting of quark and pion degrees of freedom. (We recall that this kind of double 
counting occurs, for instance, in models of hadrons based on the linear-sigma-quark- model type 
lagrangian [7,8].) This also means that we do not need to use such an ambiguous procedure 
as convoluting the pion structure functions with pion probability function (or more precisely a 
light-cone momentum distribution of the pion) inside the nucleon [9-11]. 

Several group have already attempted to calculate nucleon structure functions within the 
CQSM or the Nambu-Jona-Lasinio (NJL) sohton model. For instance, Weigel et al. investi- 
gated the polarized as well as unpolarized structure functions of the nucleon under the so-called 
"valence quark approximation" [12]. This is not an extremely bad approximation, but it is 
known to have several unpleasant features. Probably, most serious would be the violation of 
positivity condition for the unpolarized antiquark (or sea quark) distribution functions. Al- 
though such an apparent disaster does not happen for the spin dependent quark distribution 
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functions, a lesson learned from the above observation is that a reliable prediction of anti- 
quark distributions would not be obtainable unless incorporating effects of Dirac sea quarks or 
equivalently vacuum polarization effects. 

More consistent calculation including vacuum polarization effects have been performed by 
Diakonov et al. [13,14] and also by Tanikawa and Saito [15] with different regularization 
schemes, but by confining to the isosinglet unpolarized as well as isovector longitudinally 
polarized distribution functions, which have values at the leading order of l/N^. expansion (or 
at the 0th order of the expansion in the collective angular velocity Q of the hedgehog soliton) . 
Unfortunately, an abundance of interesting physics like the physics of "nucleon spin contents" 
is contained in the next order of 1/Nc expansion [5]. This is easily understood because the 
inclusion of 0{fl^) terms is the minimum condition for the collective quantization treatment 
of hedgehog solitons to hold. Otherwise, the nucleon cannot have correct quantum numbers 
[4-6]. 

We have recently reported the first calculation of the 0{fl^) contributions to the isovector 
unpolarized quark distribution function related to the physics of Gottfried sum [16] with full 
inclusion of the vacuum polarization effects [17]. It was shown that the model can explain the 
excess of the d sea over the u sea in the proton very naturally [17-19]. However, some of the 
treatments there were criticized in a recent paper by Pobylitsa ct al. [20]. In the process of 
obtaining theoretical quark distribution functions, we need to evaluate nucleon matrix elements 
of quark bilinear operators which are nonlocal in time. Their criticism is that the calculation 
in [17] docs not treat this nonlocality in time to the full extent. 

Now the purpose of the present paper is to carry out a systematic calculation of all the 
twist-2 spin dependent quark distribution functions of the nucleon as consistently as possible. 
We evaluate both of the 0{fl^) and 0{fl^) contributions with full inclusion of the vacuum 
polarization effects. The above-mentioned nonlocahty effects are also carefully taken into 
account. We believe that these unique features of our theoretical analysis would give new and 
important information on the nonperturbative aspect of the spin dependent quark distribution 
functions including the antiquark distributions as well. 

The plan of the paper is as follows. For completeness, we give in sect. 2 a precise definition 
of twist-2 quark distribution functions which we shall investigate in the present paper. How to 
evaluate these quark distribution functions within the framework of the CQSM is explained in 
sect. 3. Sect. 4 is devoted to the discussion of the numerical results. We then summarize what 
we have found in sect. 5. 
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2 Definition of quark distribution functions 



Most theoretical analyses of quark distribution functions of the nucleon are based on a 
field-theoretical formulation given by Collins and Soper [21]. As a natural extension, Jaffe and 
Ji recently carried out a systematic classification of quark distribution functions by including 
chiral-odd distribution functions which do not appear in the formulas of deep inelastic scatter- 
ing cross sections [22]. According to them, there are nine independent distribution functions, 
from twist 2 to twist 4. Here we are interested in the twist-2 distribution functions, which 
are known to have simple parton model interpretation. There are three twist-2 distribution 
functions, the spin independent (or averaged) distribution fi{x), the longitudinally polarized 
distribution gi{x), and what is called the transversity distribution hi{x). Following the nota- 
tion of [22], they are represented as 

Mx) = -^J^e^'^<PS\^ljim^{Xn)\PS>, (1) 

g^(x) = -^J^e''^<PS,\^i(OhMXn)\PS,>, (2) 

h^(x) = -l-J^e^^^<PS^\^l;i(0)^^^,M^n)\PS^>, (3) 

where and n'^ are two light-like (null) vectors, having the properties, 

p- = 0, n+ = 0, 71^ = 0, p-n^l. (4) 

Without loss of generality, one can choose a frame in which the four-momentum P^ of the 
initial nucleon and the four-momentum transfer from a lepton to a nucleon have the third 
and the time components only. In this frame, p'^ and n'^ take the form : 

p^ = -^(l,0,0,l), n'^=^ (1,0, 0,-1), (5) 

while P^ and are represented as 

P" = p^ + —n'', (6) 



- <Ju' + M%Q' ) p" + i (i' + + MJQ^ ) n^ (7) 



with 1/ — P-q and Q'^ = —q^. In the above definition of the twist-2 quark distribution functions, 
^+ is a component of the quark field -0 defined through the decomposition 

= (P+ + P_)7/> = 7/.+ + 7/._, (8) 

by the projection operators P± — §7^7^ with 7^ = 73 ^'^^ According to the authors of 

[22], -0+ is called the "good" component of -0, since it describes an independent propagating 
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degrees of freedom in the light-cone quantization scheme [23]. On the other hand, ip- is 
called the "bad" component, since it can be interpreted as quark-gluon composites. It is 
important to recognize that only the good component of ip appears in the definition of twist-2 
quark distribution functions in conformity with the fact that they have simple parton model 
interpretation. In the actual model calculation of these distribution functions, it is more 
convenient to rewrite the above expressions with use of the identities : 

Pi - - ^(I + tV), (9) 

^+75^+ = ^(1 + 7V)75, (10) 

P+^±^,P+ = ^(1 + 7V)7±75- (11) 

Since the distribution functions are in principle frame-independent, it is also convenient to go 
to the nucleon rest frame, in which one can set V = Mn/\/2. Now using the change of variable 

as 

^n^^X^il,0,0,-l) = z^, (12) 



we obtain 



Noting that 



dX e'^ • • • = Mjv / dzQ e^^^^ (14) 

-OO J — OO 

we are then led to the following expressions : 

h{x) = j dz'^e'^^'''^ 

X <P = 0,5|^t(o)(i + /^3)^(;,)|P = 0,5> |,3=_,„,,,=o, (15) 

47r J 

X < P = 0, 5, I V'^(0)(1 + 7V)75V'(^) I P = 0, 5, > |,3=_,„, ,^=0, (16) 
hi{x) = ^ [ rf^^e^^^'^"" 



47r 

X <P = 0,S^\ V^t(o)(l + 7V)7±75V'(^) \P = Q,S^> U3=_,„, ,,=0- (17) 

What is left for us now is to evaluate nucleon matrix elements of quark bilinear operators 
containing two space-time coordinates with light-cone distance. How to evaluate these matrix 
elements of bilocal quark operators will be explained in the next section. 
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3 Theory of quark distribution functions 

As shown in the previous section, the quark distribution functions of our present interest 
can generally be represented in the form : 

1 /■oo _ 

q{x) = ^ dzo e-^-^« {N{P = 0) | V'(O) rV^(z) | N{P = 0)) |.3=-.o,..=o ■ (18) 

In the present study, we confine to spin-dependent distribution functions, so that we are to 
take 

Oa = (1 + 7V)75, r3(l + 7V)75, (19) 

respectively for the isoscalar and isovector parts of the longitudinally polarized distribution 
functions, whereas 

Oa = (1 + 7V)7±75, Ta (1 + 7V) 7±75, (20) 

for the isoscalar and isovector parts of the transversity distributions. We recall here the fact 
that, extending the definition of distribution function q{x) to interval —l<x<l, the relevant 
antiquark distributions are given as [14], 

Au{x) + Ad{x) = Au{-x) + Ad{-x) (0 < x < 1), (21) 
Au{x)-Ad{x) = Au{-x) - Ad{-x) (0 < x < 1), (22) 

for the longitudinally polarized distributions, while 

Su{x) + Sd{x) = - [Su{-x) + Sd{-x)] (0 < x < 1), (23) 
5u(x) - 5d(x) = - [5u(-x) - 5d(-x)] (0 < x < 1), (24) 

for the transversity distributions [22]. As explained in the previous paper [17], the basis of 
our analysis is the following path integral representation of a matrix element of an arbitrary 
(bifocal) quark bilinear operator between the nucleon states with definite momenta : 



{N{P)\ij\0)Oij{z)\N{P)) = ^ J d^x d^y e-'^-^ e'^ y j VU 
X jVi/jVi/j^ m'^,^) ^P\0)OxIj{z) 4(-|,2/) exp[z J d^xi^{i^ - MC/^«)V], (25) 
where 

C = ^ - MW{x))i) , (26) 

with U^^{x) = exp[i75T • 7r{x)/ f-;^] being the basic lagrangian of the CQSM, and 

Mx) = ^^""••""^ry;3;/T3^^„,;,(x)---t^«,^;,^(a;), (27) 



is a composite operator carrying the quantum numbers JJ^^TT^ (spin, isospin) of the nu- 
cleon, where ccj is the color index, while ^Yj^^tt^^ is a symmetric matrix in spin-flavor in- 
dices fi. By starting with a stationary pion field configuration of hedgehog shape Uq^{x) — 
exp [ i75T • rF{r) ] , the path integral over the pion fields U can be done in a saddle point 
approximation. Next, we consider two important fiuctuations around the static configuration, 
i.e. the translational and rotational zero-modes. To treat the translational zero-modes, we 
use an approximate momentum projection procedure of the nucleon state, which amounts to 
integrating over all shift R of the sohton center-of-mass coordinates [14] : 

{N{P)\ijj\0)Oilj{z)\N{P)) — > J d^R {N{P)\il:\0,-R)Oilj{zo,z- R)\N{P)) . (28) 

The rotational zero-modes can be treated by introducing a rotating meson field of the form : 

U^^{x,t) = A{t)U^'{x)A\t) , (29) 

where A{t) is a time-dependent SU{2) matrix in the isospin space. Note first the identity 

^- MA{t)U^'{x)A\t))ip = ip\{idt- H -n)iPA (30) 

with 

V'A = A\t)ilj , H = + M(3U^'{x) , O = -iA\t)A{t) . (31) 

Here if is a static Dirac Hamiltonian with the background pion fields Uq''{x), playing the role 
of a mean field for quarks, while Q — \ is the SU(2)-valued angular velocity matrix later 
to be quantized as fi^ ^ Jq// with I the moment of inertia of the soliton and Ja the angular 
momentum operator [4-6]. We then introduce a change of quark field variables ip ^ ipA, which 
amounts to getting on a body-fixed rotating frame. Denoting ipA anew as ip for notational 
simplicity, the nucleon matrix element (25) can then be written as 

{N{P)\iP^{0)Oilj{z)\N{P)) 
^ 1 r{/} r^^}* J d'x d'y e-^P ^ e^P y j d'R 

r f Nc rp rp 

X / VAVipViP^ exp[i / d'^x ip^idt- H [] [A{-) iPf.{-,x)] 

J J ■ -I 2 2 

1=1 

Nc rp rp 

X i^\Q-R) A^{Q)OA{zo) i^{z,,z-R) J{ [i^i{--,y) A^{--)] ■ (32) 

Now performing the path integral over the quark fields, we obtain 

{N{P)\i)\Q)Oi}{z)\N{P)) 

1 fin fM^ j d^'x d^y e-'P-"" e'^V J d'R J VA 



Z 
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X { /. (l- ^ I i9,-H-n I ■ (A\o)o.A(^«)U ■ ^ - « I is,-H-n I - !• 

- ( - - « I a,-'„.n I °- -«>-4'(0)O,^(.„) ) ;,(|, X I I - |, } 

Nc rp ■ rp 

^ n [/.(^.^l ,g 2'^^^^] • exp[Splog(ia,-//-0)] , (33) 

with f = r{^> )] ' etc. Here Tr is to be taken over spin-flavor indices. Assuming a 
slow rotation of the hedgehog soliton, we can make use of an expansion in Q. (Since Q is 
known to be an 0{1/Nc) quantity, this perturbative expansion in Q can also be taken as a 
1/Nc expansion. For an effective action, this gives 

Sp\og{idt-H -n) = Sp\og{idt-H) + i^I J flldt. (34) 

The second term here is essentially the action of a rigid rotor, which plays the role of the 
evolution operator in the space of collective coordinates. We also use the expansion of the 
single quark propagator as 

f T i i 

- J dz'odh' f^{-,x\-^—j^\z'o,z')a -in^piz'o)- p{z'o,z'\-^—j^\0,-R)^ 

+ ... . (35) 

An important suggestion made in a recent paper by Pobilytsa et al. [20] is that one must 
also take account of the nonlocality (in time) of the operator A^{0)OaA{zo)- Expanding this 
operator around or one respectively obtains 

A^{0)OaA{zo) = At(0)O„A(0) + ^o^^(0)0,i(0) + (36) 
or A\0)OaA{zo) = A\zo)OaA{zo) - zoA\zo)OaA{zo) + (37) 

Since both choices are known to lead to the same answer [20], it is convenient to use a sym- 
metrized form in the following manipulation. This amounts to performing the following re- 
placement : 

A\0)OaA{zo) A^OaA + ^zoiA^OaAA^A-A^AA^OaA), 

^ da + iZo^{n,da}, (38) 

in the process of collective quantization of the rotational motion. Here we have introduced the 
notation 

da = A^OaA, (39) 
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for saving space. Eq.(38) means that the nonlocahty of the operator A^{0)OaA{zQ) causes a 
rotational correction proportional to the collective angular velocity il. After taking all these 
into account, we are then led to a perturbative series in which is also regarded as a 1/Nc 
expansion : 

(7V(P)|V;t(0)O„V(^)|iV(P)) 
= {N{P)\^^;\0)Oa^{z)\N{P)f + {N{P)\^^;^ {0)Oa^{z)\N{P)f + (40) 

where 

{N{P)\il^\0)Oil>{z)\N{P)f° 



X 



Z ^ ^ y . . y o o J u... J y^a,,S 

-s{zo,z-R\j^\0,-R),.f,{-,x\j^\ --,y),. 



N, 

X 

J=2 



Nc rp ■ rp J 

n t/^-^2'^'ia^' " 2'^^^^-^ ■ exp[Splog(za,-i/) + z-yf)^dt] , (41) 



and 



{N{P)\i^\Q)0^l:{z)\N{P)f' 



X |y dS' dz', iQa/siz',) (A^(0)naA(zo))^s 



X 



-fA^,x\-^\-^,y),,.s{zo,z-R\^^\ziz\.,{z',,z'\-^ 

+ izo-{n,da}^s 



X 



X 11 " 2'^^*^ ■ exp[Splog(^a,-if)+.-|(]^cit]. (42) 

j=2 * 
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Let us first discuss the leading 0(f2°) term. As usual [4,5], we introduce the eigenstates \m) 
and the associated eigenenergies Ej^ of the static Dirac hamiltonian H, satisfying 

H\m) = E^\m). (43) 

This enables us to write down a spectral representation of the single quark Green function as 
follows : 

- 0{t'-t) Yl e-''''-'^'-''\{x\m){m\x')p . (44) 

m<0 

Using this equation together with the relation 

{z-R\ = {-R\e'P-^, (45) 

we can perform the integration over R in (41). The resultant expression is then put into (18) 
to carry out the integration over zq. We then arrive at a formula, which provides us with 
a theoretical basis for evaluating the zeroth order contributions in Q to quark distribution 
functions of the nucleon : 

g(x;QO) = I ^%[U 0^'\U] dU, (46) 

where 

= (-1)""" ■^-■^.^.K-) ■ (*^) 

are wave functions, describing the collective rotational motion of the hedgehog soliton, while 
0^'\U] = Mr,^{Y-E)i^\OaS{xM^-Er,-p')\n). (48) 

^ n<0 n>0 

Using the identity 

( E + E ) (^1 OaSixMj, - K - p') \n) = 0, (49) 

n<0 n>0 

eq.(48) can be expressed in either of the following two forms : 

0(0) [Ca] = M^N^Y {n\ OaS{xM^ - E^ - p^) \n) 

n<0 

= - Mjv A^c E (^1 OaS{xMN -En- p') \n) , (50) 

n>0 

i.e., as a sum over the occupied states or as a sum over the nonoccupied states. As was 
emphasized in [14], it is better to use the first form for x > 0, whereas the second form for 
x < 0, for the purpose of numerical calculation. 
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Next we turn to the 0{fl^) contribution. In writing down (42), we have retained the 
time arguments 0,^0 and z'q in ^4^, A and since we have to pay attention to the time 
order of these collective space operators, which do not generally commute after collective 
quantization of the rotational zero-energy modes. In the previous paper [17], motivated by the 
physical picture that the time-scale of deep inelastic-scattering processes is much shorter than 
that of collective rotational motion of the soliton, we dropped special time-order diagrams 
in which the Coriolis coupling Q between the collective rotational motion and the intrinsic 
quark motion operates in the time interval between zq and 0. However, this procedure was 
criticized by Pobylitsa et al. in a recent paper [20]. According to the them, there is little 
reason to assume approximate degeneracy of and Zq in A^{0)OaA{zo), since the deep-inelastic 
scattering processes are not necessarily short distance phenomena. Taking this nonlocality in 
time arguments more seriously, one should retain all the possible time-order diagrams. In 
doing so, we must pay attention to the time order of collective space operators A and fi. By 
ordering these operators according to their time orders, we are led to the replacement : 

Q^f,{z',) {A^{0)OaA{zo)\, 

— > [9{zq,0,zo) + 9{zQ,zo,0)]Qapdy5 + [0{0,zo,Zq) +9{zo,0,ZQ)]d^sQap 
+ e{0,z'o,zo){Oa)^,s,A\^,Q^f3As>s + e{zo,ziO){Oa)^,s,As>sn^/3A\y. (51) 

Here the third and the fourth terms are new ones discarded in the treatment of [17]. In order 
to handle these somewhat peculiar terms, we first recall the rule of collective quantization : 

Q = ^ QaTa " ^ JaTa, (52) 

where Ja is the total angular momentum operator satisfying the commutation relations (CR) 
as follows : 

[Ja, Jb] — i ^abc Jc, (53) 

[Ja,A] = ^Ara, (54) 

[Ja,A^] = -l^aAl (55) 

Using these CR, one can show that 

{Oa)^>S' A\y ^a(3 As'S = — (Tc)^^ {Oa)^>S' ^77' '^'^^^'^ 

= ^(^c)«^[^(^^0„^re),,+ (^to„^),,Je], (56) 
where we have used (54). Similarly, by using (55), one may obtain an alternative expression 

{Oa\., A\y n^fs As>s = Yi [ \ i^cA^O^^U + Jc {A^O,A)J. (57) 

11 



In the following manipulation, we find it convenient to take an average of these two expressions 
as 

+ ^j{rc)^p[{A^OA)^,J, + UA^O,A\,]. (58) 

Now we must treat two cases separately. The first is the case in which the operator Oa contains 
an isospin factor Ta as 

Oa = TaO. (59) 

In this case, using the relation A^OaA — DabTbO, we can rewrite as 

{A^OaAr,)^s + {rcA^OaA)^g = D^b {{nrc + rcn)0)^^ = 2 D^c (O)^,. (60) 
On the other hand, if Oa contains no isospin factor as 

Oa = O, (61) 

we obtain 

{A^OaAT,)^,+ {rcA^OaA)^, = 2 (tcO)^,. (62) 
Unifying the two cases, we can then write as 

{Oa)ys' A\y n^pAs^s = ^ I 1^";^^"^ I + I {^a,, (A^OaA)^,}^. (63) 

A similar manipulation for the fourth term in (51) leads to 



iOa)ys'As'Sn^pA\^. = - ^ {rX, I ^^^^"^ I + \ {^ap, (^^0„A),Jh 



(64) 



Retaining all these possible time order diagrams, the O(Q^) contribution to the distribution 
function now becomes 



{N{P)\ilj^z)Oilj{0)\N{P)f 
Z 



f^^> f^^>' I d^x d'y e-*^-^ e'Py J d^R J VA 



X {i I d^z'dz'o 



X [^(4, 0, 2:0) + ^(4, zo, 0)] QapO^s + [0(0, zo, z'^) + e{zo, 0, 4)] O^s^a. 

+ 0{O, 4, ^0) [ li^ap, 0,s}^ + ^(re)„, I ^^^^"'^ 
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+ e(zo, Z'o, 0) f 0^S}+ - Jj{'^c)ap | 



X 



T i 



idt-H 
i 

idt-H 



idf-H 



0, —R)j ■ s{zo, z — R \ - 



0, —R)f ■ s{zq, z — R\- 



idf-H 



Zqi Z )a ■ Z I . 



idt-H 

i 

idt~^^ ' 



X 



idt-H 



-s{zq,z-R\- 



T 

2^y)9i ■ s{zo,z- R\ 



0, —R)^ ■ si^o, z — R\ 



Nc rp 

X n [/.(^'^ 

j=2 ^ 



idt-H 

i 



idt-H 



idt-H 

i 



Z )a ' /3(^0' ^ 



idt-H 



T 

0, -R). 



idt-H 



T 

2 ' Vlgi 

T 
2" 



idt-H 



2 ' 



eM^V^og{idt- H) + i^- j nldt] , 



(65) 



After stating all the delicacies inherent in the structure function problem, we can now proceed 
in the same way as [17] and [24]. Using the spectral representation of the single quark Green 
function (44) together with the relation (45), we can perform the integration over R,z', and 
z'q. The resultant expression is then put into (18) to carry out the integration over zq. We then 
arrive at a formula, which gives a theoretical basis for evaluating the O(Jl^) contributions to 
quark distribution functions of the nucleon : 



where 



with 



?(^;^^') = j o^'\u] ^%Ma] dU 



(66) 
(67) 



M. 



N 



M, 



4 



N 



E - E 



E - E 



E 


— En 


\n) 


+ (^1 




1 


E 


— En 


\\n) 


+ {n 




1 


Era 


- En 



X [ (n I Oa{5n - 5m)\m) {m\Vt\n) + {n\Vt\m) {m\Oa (5n - ^m) I ], (69) 



X (n I Tc I m) (m 
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DacO 



(70) 



= ^^iE-i:){n\{OaM}Sn\n). (71) 

^ n<0 n>0 

In the above equations, we have used the notation 

5m = 5(xMN-Em-p^), and 5n = 5(xMn - E^, - p^). (72) 

for saving space. Here O^"* is the contribution from the diagram in which z'q is later (or 
earher) than both of or zq. As was emphasized in [17], this term contains transitions between 
the occupied and nonoccupied single quark levels so that it is not in conflict with the Pauli 
principle. On the other hand, O^"* and O^? are the contributions from diagrams in which 
2;q lies between and zq. Although these terms appear to contain Pauli-violating transitions 
between the occupied levels themselves or the nonoccupied ones, wc take here the viewpoint 
advocated in [20] that there is no compulsory reason to drop them since we are here dealing 
with operators which are non-local in time. Finally, Oq^ is the 0(f2^) contribution resulting 
from the nonlocality of the operator A'' {O)OaA(zQ), i.e. the second term of (38). In deriving 
Oq\ use has been made of the identity, 

^ dzoi zo (^^--^r.-P^) ^0 ^ _E^_ _ (73) 

27r J-oo Mn ox 

As will become clear shortly, it is convenient to treat O^^^ and in a combined way. To 
see it, first note that, after a simple change of summation indices, can be rewritten as 

Oi'^ = Mn^{ E pip [{n\daSn\m){m\n\n) + {n\n\m){m\djn\n)] 

^ m>0,n<0 

~ E F [(^|Qa^n|^)(^|^|^) + (^|^|^)(^|Qa^n|"^)]}- (74) 

m<0,n>0 -^"^ 

From now on, wc treat the two cases separately. First, assume that the relevant operator Oa 
contains an isospin factor in such a form as Oa = TgO. In this case, in view of the relations 
Oa = A'^OaA = DabTbO and f2 = ^Jc^c, we must carefully treat the noncommutativity of 
the two collective space operators Dab and Jc- By keeping the order of Dab and Jc, O^'* can 
generally be divided into two pieces [24] as 

= M^^^{Dab,JcU 



X 



{ E p \ p [{n\TbOdn\m) {m\Tc\n) + {n\Tc\m) {m\TbOS„\n)] 

m>0,n<0 ^'^ 

E F ^-F [{'^\'^b05n\rn){m\Tc\n) + (n|rc|m)(m|rft05„|n)]} 



E - E 

m<0,n>0 -^"^ 

N 1 

+ M^^-[Dab,J,^ 
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X 



{ i: 



E — E 

m>0,n<0 -^'^ 



[{n\n05n\m){m\Tc\n) - {n\Tc\m){m\n05n\n)] 



^ — [(n|r605„|m)(m|rc|n) - (n|rc|m)(m|rbC'5„|n)]}, (75) 



£J — E 

m<0,n>0 -^"^ 



which contains symmetric and antisymmetric combinations of the two collective space opera- 
tors Dab and Jc- On the other hand, it can be easily verified that O^q term contains symmetric 
combination only : 



O 



(1) 



N 1 



X 



[{n\Ti,05n\rn){m\Tc\n) + {n\Tc\rn) {m\Ti,05n\n)] 



l ^ E - E 

m>0,n<0 -^"1 

~ E Tr^-^r [(^k60(5„|m)(m|rc|n) + (n|rc|m)(m|T605„|n)]}. 

m<0,n>0 -^n 

Combining O^"* and terms, we then obtain for the isovector case 



where 



O 



(1) 



Mn^\ {Dab. Jc}+ I E 



- E + E - E 

m>0,n<0 m<0,n>0 m<0,n<0 m>0,n>0; 



(76) 



(77) 



(78) 



X 



E — E 



{n\TbOSn\m) {mlTcln) + (n|Tc|m)(m|Tb05„|n) 



O 



(1) 

[A,B] 



iV. 1 



MN^^[Dab,Jc] I E + E 



X 



m>0,n<0 m<0,n>Oy 

{n\Tf,06n\m) {m\Tc\n) — {n\Tc\m) {m\Tb06n\n) . (79) 



Em — E,_ 

The situation is much simpler for isoscalar operators Oa — O. Since Oa — A^OaA — 
A^OA — O, we have only to replace both of Dab a-nd Tb by 1 in the above manipulation, 
thereby leading to 



O 



(1) 

{A,B} 



- Mn^jA E - E + E - E 

\m>0,n<0 m<0,n>0 m<0,n<0 m>0,n>oy 

X ^ {n\05n\rn) {m\Tc\n) + {n\Tc\m) {m\05n\n) , (80) 



— £■ 



[A,B] 



0. 



One notices that only the symmetric combination of the matrix elements survives for this 
isoscalar case. This should be contrasted to the isovector case in which either of the sym- 
metric part or the antisymmetric part survives, depending on the symmetry property of the 
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relevant single quark matrix elements appearing in (78) and (79). As we shall discuss later, the 
symmetric part contributes to the isoscalar unpolarized distribution function Au{x) + Ad{x) 
and 5u{x) + 5d{x) at the O(Q^), whereas the antisymmetric part plays an important role in the 
0{D}) term of the isovector polarized distribution functions Au{x) — Ad{x) or 5u{x) — 5d{x) 
[24]. 

Now we shall investigate the case of our interest in more detail for obtaining explicit 
formulas, which can be used for numerical calculation of polarized distribution functions of 
the nucleon. 



3.1 Au{x) + Ad{x) 

The relevant operator in this case is 

da = At(l + 7V)75^ = (1 + 7V)7'- (82) 

Since the 0{Q^) contribution to Au{x) + Ad{x) vanishes due to the hedgehog symmetry, the 
leading contribution to this distribution function arises from the 0{Q^) terms. Due to the sym- 
metry property of the relevant single quark matrix elements, only the symmetric combination 
of O^^^ + O^'' survives. The total 0{Q^) term therefore consists of three pieces, O^^J^^, Ob' 
and Oq^ . Using the general formulas obtained so far, the contributions of these three terms 
to Au{x) + Ad{x) are given as 



[Au{x) + Ad{x)]fl^y = {Js)^,-M^'-f{ E - E 



4/ ^ ^ E — E 

m=all,n<0 m=all,n>0 " 

X [(n|(l + 7%^)75(5„|m)(m|r3|n) + (n|T3|m)(m|(l + 7V)75(^n|^)], (83) 

1(1) _ /i\ . Ayr,. ^ 



[Au{x) + Ad{x)]^^] = {l),,-MN^i E - E ) 



m<0,n<0 m>,ra>0 



X {n\Tc\m){m\Tc{l + 7V)75(5n - Sm)\n), (84) 
[Au{x) + Ad(x)]« = (Ja),, -^M^^ij:-^) {nUl + ^VHW- (85) 

^■^ n<0 n>0 



In the above equations, (C')p-f denotes a matrix element of a collective space operator O with 
respect to the proton in the spin up state along the z-axis, i.e. 

{O)^^ = f^[^hU]0^[^hUdU = {p,S, = l/2\0\p,S, = l/2). (86) 

J 2 2 2 2 

In deriving (83), we have used the relation 

= ({(l + 7V)75,^JcTe}^)^^ = (J3),T -^3(1 + 7%') 75. (87) 
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One may notice that the collective space operator contained in the term [Au{x) + Ad(x)]^^} 
is 1 and it is different from J3 contained in other two terms. The appearance of this term 
seems to be inconsistent, since it does not change sign in contrast to the other two terms when 
the direction of the proton spin is reversed. Fortunately, it can be shown that this potentially 
dangerous term vanishes identically due to the symmetry of the double sum of the single quark 
matrix element : 

[Au{x) + Adix)]^^} ^ 0. (88) 

We are then left with the two terms, i.e. [Au{x) + Ad{x)]^^^^y and [Au{x) + Ad{x)]^(^\ which 
both have required state dependence. For the purpose of numerical calculation, it is convenient 
to rewrite the above two terms slightly further. Using the argument given in appendix of [15], 
we can prove the identity : 

( E + E 



E - E 

m=all,n<0 m=all,n>0 

X [(n|(l + 7%2)75(5„|m)(m|r3|n) + (n|T3|m)(m|(l + 7%2)75(5„|n)] = 0. (89) 
As a consequence, [Au{x) + Ad{x)]^^^l can be expressed either of the following two forms : 
[Au{x) + Ad{x)]fl^^ 

N 1 

= i-^^lr^N-r^ E 



m=all,n<0 -^"^ 

X [(n|(l + 7%^)755n|m)(m|T3|n) + {n\T3\m){m\{l + -f^-f^)-f55n\n)] 



pT ^ at E - E 

m=all,n>0 

X [(n|(l + 7%^)755„|m)(m|T3|n) + (n|T3|m)(m|(l + 7%^)755„|n)]. (90) 

As advocated in [20], it is convenient to use the first expression given as a sum over the 
occupied states for the numerical calculation of distribution functions in the region x > 0, 
while to use the second one given as sum over the non-occupied states when x < 0, since one 
can thus avoid vacuum subtraction, i.e. subtraction of the corresponding sums over vacuum 
levels (with U = 1). Following [20], we also separate the Em = En contribution from the above 
sum over the single quark levels. This can be done by noting the identities, 

E F i F (n|(l + 7V)75(^n|m)(m|T3|n) + (n|r3|m)(m|(l + 7V)75(^n|n)] 

m<0,n<0 -^"^ 

= ^ E F j_ F [(^l(l + 7V)75((^n -(^m)|m)(m|T3|n) 

^ m<0,n<0 

+ (n|r3|m)(m|(l + 7V)75(<^n-UI^)], (91) 
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and 

^.^ 5{xMm -En- p") - 6{xMn -Em- p") 

Em — >En Ejn — En 

= S'ixMN-En-p") = -^^S{xMj,-En-p'). (92) 

Mn ax 

Prom (90), we can then readily obtain 

(1) 



[Au{x) + Ad{x)]\l^^ 



pT ^ at ^ E - E 

m=all,n<0 "^^ 



X 



(n|(l + j"j'^y/5Sn\m){m\T3\n) + (n|T3|m)(m|(l + 7"7^)75'5n|«) 



^ ^ ''P^ dx 4/ i-^^ - En 

m<0,n<0 

{Em=En) 



X 



(n|(l + 7°7=^)75(5„ - 5„)|m)(m|T3|n) + (n|T3|m)(m|(l + 7°7=^)75(5„ - 5m)\n) , (93) 



and a corresponding expression given as sums over non-occupied levels. The remaining term 
can similarly be expressed in either of the two equivalent forms as 

[Au{x) + Ad(x)]« = ( Ja),^ -Ijfi: (^1^3(1 + 7V)755n|n) 

"•^ n<0 

= - ^■^3)pT • 77 E (^1^3(1 + 7V)755„|n). (94) 

n>0 

Inserting the complete set of single quark states into the first expression and separating the 
Ejn 7^ En and E^ — En terms in this sum, we obtain 

[A«(a;) + Arf(a;)]g) = {Js),^ ^JJ E {n\rs\m) {m\{l + ^VHSn\n) 

•m=all,n<0 
(Em^En) 

+ E (n|r3|m)(m|(l + 7V)75<^n|n), (95) 

m<0,n<0 

(-£'m=-£'n) 

and a corresponding expression given as sums over non-occupied states. Just as argued in [15] 
for the case of the unpolarized distribution function u{x) — d{x), Em = En contribution in 
the double sums in [Au{x) + Ad{x)]^^l^^y and [Au{x) + Ad{x)]^(p precisely cancel each other. 
After regrouping the terms in such a way that this cancellation occurs at the level of analytical 
expressions, the 0{fl^) contribution to the distribution function Au{x) + Ad{x) can finally be 
written in the following form : 

[Au{x) + Ad{x)f^ = [Au{x) + Ad{x)]fl^y + [Au{x) + Ad{x)]^^} (96) 
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where 

1(1) 



[Au{x) + Ad{x)][l^y 
= (^3>,T-^iv§ E E -E (^l(l + 7V)75<^n|m)Hr3|n), (97) 

m=all,n<0 

= - ( Ja)^^ • Mjv ^ E p; i (n|(l + 7V)75^n|m)(m|r3|n), (98) 

m=aU,n>0 " 



and 



[A«(x) + Ad(x)]g) = ^'^'^prj^JJ E (r^|T3|m)(m|(l + 7V)755n|n), (99) 

m=aii,n<0 

= -(■/3)pT-^§ E (n|r3|m)(m|(l + 7V)755n|n). (100) 

m=aii,n>0 

(-E-m ) 

These expression will be used in the numerical calculation. 



3.2 Au{x) - Ad{x) 

The relevant operator for the isovector longitudinally polarized distribution function is 

Oa=3 = AtT3(l+7V)75^ = i?36rb(l + 7V)75- (101) 

The main contribution to this distribution function comes from the 0th order term in Q. A 
simple manipulation gives 

[Au{x) - Ad{x)f^ = (Das)^^ -MnN^YI (MMl + iVHSnln) 

n<0 

= -{Dss)p^-Mr,N, E (^1^-3(1 + 7V)75^n|n). (102) 

n>0 

The 0{Q^) contribution to Au{x) — Ad{x) is much far comphcated. It generally consists of 4 
terms, 0'^\ Og? and O^c^ . As was already mentioned, the symmetric part of the sum of 
O^^ and vanishes for this particular operator, owing to the symmetry of the single quark 
matrix elements. Using the familiar commutation relation 

[Jc,Dsh\ = iecbeD3e, (103) 

the antisymmetric part of O^^^ + becomes 
[A«(x)-A(i(a:)]«^j 
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Nc 1 



E — E 

m>0,n<0 -^"^ 



X [{n\Tc\m) {m\n{l + 7V)75((^n + 5^)\n) + (n|rfe(l + 7V)75('5n + 5,n)\rn) {m\Tc\n)] 



^ ^ m>0,n<0 -^"1 



X [(n|r+i|m)(n|r+i(75 + E3) 



5n + 



m) - (n|r_i|m)(n|r_i(75 + E3) 



m)], (104) 



with the standard definition t± — ^(ti ± iT-^j\pl and E3 = 7"7'^75. Next, from (70) with the 
case of isovector operator, we find that 



1 



E - E 

km<0,n<0 m>0,n>0; 



X 



{n\T3\m){m\{l + - 5m)\n). (105) 



E - E 



One should notice that the state dependence of this somewhat pecuhar contribution is nothing 
different from that of the main term, which imphes that there is no reason for this term 
to vanish. In fact, the single quark matrix element appearing in the above double sum is 
essentially the same as that appearing in the expression for [AM(a;) + Ad{x)]^^\ 

The last but potentially important contribution comes from the nonlocality correction term 
Oq\ First note that 



1 



= ^ [ 2 ^^3^' -^'^^ + iVHtc + Tcnil + 7%')75 
+ ^ [^36, Jc] [nil + 7°7')75rc - + 7°7')75 

^ ({i?3c, Jc} + (1 + 7°7')75 + i ebce [D3b, Jc] Te (1 + 7%')75) • (106) 
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The first term of the above quation does not contribute, since 

J2 {n\{l + iVHSnln) = E (n|(l + ^VHK\n) = 0. 

n<0 n>0 

Simplifying the second term by using the CR (103), we finally obtain 



(107) 



[Au{x) - Ad(x)]g' = - (D33)jy • ^ § I E - E 1 (^1^3(1 + ^Vh%\n). (108) 



in<0 n>0; 



For the same reason as before, it is convenient to consider these term in a combined way. 
To this end, we first rewrite (108) by inserting a complete set of single quark states and by 
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separating the — En contributions from the resultant double sum. The result can be 
expressed in two alternative forms as 



[Au{x) - Ad{x)f^ 



c 



(D ) ^ 



\ 



E + E 



(n|T3|m)(m|(l + 7 7 )^-,5n\n) 



{D33) 



d_N, 
' dx 41 



m=all,n<0 m<0,n<0 

\ 



E + E 

m=all,n>0 m>0,n>0 
\ Em^E-n E^= 



{n\n\m) (m| (1 + 7V)75<^n|n) . (109) 



To rewrite the B' term, we first separate E^ — En contributions in the double sum of (105) as 

[Au{x) - Ad{x)]^^} 



I 



E - E 



m<0,n<0 m>0,n>0 

\{Em^E„) (Em+E-a)) 



X 



E - E 



(n|T3|m)(m|(l + 7 7 )75(^„|n), 



+ (-/^33),--- 



\ 



E - E 



(n|r3|m)(m|(l + 7%3)^5'5„|n). (110) 



TO<0,n<0 m>0,n>0 

\{Em=ErC) {Em=En) / 



Next, we notice the identity 

= E 

1 d 



{n\T2\m){m\{l + 7°7^)755„|n 



m=all,n=all 

(Em^En) 



^ 2M^di ^ (n|T3|m)(m|(l + 7V)75<^nk). 



m=all,n=all 
{Em=E„) 



Using this identity, (110) can be rewritten in either of the following two forms : 

[Au{x) - Ad{x)]^s' 



I 



N 



\ 



m<0,n<0 m>0,n<0 m<0,n>0 

\ (EmJ^En) I 

X p ^ p {n\r3\m){m\{l + -f^-f^)-f55n\n), 



(111) 
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d N, 



m<0,n<0 

{Em=En) 



N 



2 E + E + E 

m>0,n>0 m>0,n<0 m<0,n>0 

X———— (n|T3|m)(m|(l + 7%''^)755„|n), 



d N, 



-(^33)pt-^^ E HT3|m)(m|(l + 7V)75<^n|ri), 



(112) 



m>0,n>0 

{Em=En) 



Comparing (109) and (112), one notices that the Em = En pieces in the double sums cancels 
precisely between B' and C terms. (This is true for both of the occupied and nonoccupied 
expressions.) After some manipulation by taking care of this cancellation, the sum of these 
two terms can finally be expressed as 



[i^u{x) - M{x)]%}^c 

{ E w\w HT3|m)(m|(l + 7V)75(^nk), 



X 



m=all,n<0 "^^ 

(Em^En) 

1 d 



2Mn dx 



(n|T3|m)(m|(l + 7V)75^n|?^), 



m=all,n<0 

{Em¥=En) 



E 

^ E - E 

m>0,n<0 -^"1 



(n I T3 1 m) (m I ( 1 + 7 V) 75 ^^^^^-^ I } 



X 



I E 

m=all,n>0 ^"^ "^^ 



{n\T^\m){m\{l + 7°7^)755n|w), 



1 ci 



2Mn dx 

m<0,n>0 



J2 (n|T3|m)(m|(l + 7%^)755„|n), 



m=aZZ,n>0 

(Em^En) 
1 



{n\Tz\m){m\{\ + 7^7^)75 



|n)}. (113) 



For numerical calculation, we shall use the first form for a; > 0, while the second form for 
a; < 0. 



22 



3.3 6u{x) + 5d{x) 

Since the evaluation of the transversity distribution can be done in a completely parallel way as 
the longitudinally polarized distribution functions, we shall show below only the final results. 
The 0(0°) contribution to Su{x) + Sd{x) vanishes, i.e. 

[Su{x) + Sd{x)f^ = 0. (114) 

The 0{Q^) contribution consists of two pieces as 

[6u{x) + Sd{x)f^ = [Su{x) + 6d{x)]flj^y + [6u{x) + Sd{x)]'§} , (115) 

where 

[Su{x) + Sd{x)]^^lsy 

= WpS:.-^N^ -^^r^{'n\{l^l5-il'^)Sn\m){m\Ti\n), 

m=aU,n<0 " 

^ -('^^)ps^- TP ^ TP {n\{-f^-f5 - h'^)d„\m){m\Ti\n), (116) 



and 



m=all,n>0 " 



[Aii(x) + M{x)]^S 

m=all ,n<.0 

(Em^En) 

= "^"^^W'^^ 51 {n\n\m){m\{-f^-f5-i-f'^)Sn\n) (117) 



m=all,n>0 

{Em^En) 



Here {Jx)ps^ is defined by 



3.4 ^^(a:) — 6d{x) 

The 0(^1°) contribution is given by 

[6u{x) - 6d{x)f^ = {D,,)^s. -M^N.Y {n\T^{i^l, - ^7')4|n) 

n<0 

= - (i^3l)p5. -MnN^Y. («k3(y75 - i7')<^n|n). (119) 

n>0 
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The 0{il^) contribution consists of two pieces as 

[Su{x) - Sd{x)f^ = [Su{x) - Sd{x)]\2s] + IH^) - Hx)]b'+c (120) 

where 

[5u{x) - Sd{x)f2,B] 

N. 1 

°^ m>0,n<0 ^'^ 

X [(n|Tc|m)(m|Tb(7S5-i7^)(5„ + ^m)!"-) + {n\Tb{-f^-f5 - i-f'^){Sn + Sm)\m) {m\Tc\n)] 
IN 1 

^ m>0,n<0 

x[(n|r+i|m)(n|r+i (7^5-^7^) ^" ^ ^"" [m) - (n|r_i|m)(n|r_i(7^75 - i7^)^^^-t-^|m)](121) 

and 

[(5m(x) - Sd{x)]^^}^c 

N 

1 



^ { ^ F F (^k3|m)(m|(7S -n^)^n|^) 

- 2i^^ E (^k3|m)(m|(y75-^f)<^n|n) 

m=all,n<0 

E F (^k3|m)(m|(7S-n^) ^"^'^"» }, 

m>0,n<0 -^"i ^ 

TV 

= -(^3lU-^iV^ 

{ E F j_ F (^k3|m)(m|(7S -^7^)'^nk), 



m=a/i,n>0 
1 d 



- S (^k3|m)(m|(7S5 -i7^)5n|n), 

(£'m.7^£'n) 

- E Tr^Hr3|m)(m|(y75-^7^)^^|r^)}. (122) 

m<0,n>0 ^ 
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4 Numerical results and discussion 



Before showing the results of numerical calculations, we briefly discuss the parameters 
of our effective model specified by the lagrangian (26). Fixing to its physical value , i.e. 
/tt = 93 MeV, only one parameters of the model is the constituent quark mass M, which plays 
the role of the couphng constant between the pion and the effective quark fields. There is 
some argument based on the instanton picture of the QCD vacuum that the value of this 
mass parameter should not be extremely far from 350 MeV [25]. Phenomenological analyses 
of various static baryon observables based on this model prefer a slightly larger value of M 
between 350 MeV and 425 MeV [5,6]. In the present analysis, we use the value M = 375 MeV 
favored from analyses of various static observables of baryons. Actually the model contains 
ultraviolet divergences so that it must be regularized by introducing some physical cutoff. In 
the case of static nucleon observables, most frequently used regularization scheme is the one 
based on Schwinger's proper-time representation [5,6]. Unfortunately, how to generalize this 
regularization scheme in the evaluation of nucleon structure functions is an open problem. For 
evaluating quark distribution functions, Diakonov et al. then proposed to use the so-called 
Pauli-Villars regularization scheme, which they claim has several nice properties as compared 
with the energy cutoff scheme like the proper-time regularization scheme [14]. The basic idea 
of this regularization scheme is very simple. Using the derivative (gradient) expansion, one 
can evaluate the effective meson action corresponding to the original effective quark lagrangian 
(26) as 

= ^HM).^-{d,7vf + ••■ . (123) 
Here the coefficient of the pion kinetic term given by 

- '/(2^(*J31?5yi. (124) 

contains logarithmic divergence. Clearly, this divergence can be removed by introducing a 
regularized action S^^f by 

/ A/f \ ^ 

Here S^j^ denotes the effective meson action obtained from S^j by replacing the dynamical 
quark mass M with the Pauli-Villars mass Mpy. In fact, this replaces hiM) with 
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which is clearly finite. Demanding further that the pion kinetic term in S'g^^ has the correct 
normalization, one obtain 

I^A/Mog^^/J, (127) 

For M = 375 MeV, for instance, this gives Mpv — 562 MeV. Other obscrvablcs like quark 
distribution functions, which contains logarithmic divergence, can similarly be regularized as 

M ^2 
fpy. 

For the sake of consistency, a soliton solution should also be obtained in the same regularization 
scheme. The startingpoint of soliton construction is the mean field equation 



(Or - {Of-^-^)\of-^. (128) 



(V^^)7 sinF(r) = (^^75T^)^ cosF(r), (129) 

obtained under the assumption of the static hedgehog configuration 

7r{r) = UrF{r). (130) 

Here (i'lpY/^ and {i> i jsTtpY^^^ are the regularized scalar and pseudoscalar densities in the 
Pauh-Villars subtraction scheme : 

^ (M)'' - ( 

Mpv 



(in^T-i,)':"' ^ itn^T^)" - [-^) t^i.ii.T-^t'-' . (132) 

Recently, self-consistent solutions of this equation of motion has been obtained in [26] with use 
of the Kahana-Ripka basis [27]. (Essentially the same equation was solved in [28] within the 
framework of the Nambu-Jona Lasinio model with an ad hoc nonlinear constraint.) However, 
one should use this regularization scheme with some care. In fact, it is known that the single 
subtraction is not enough to get rid of linear divergences, for instance, contained in the expres- 
sion of the vacuum quark condensate [28] , which implies that < ipip >l^^ and < i/ji^y^Til) 
also contain convergences. Why could the authors of refs. [26,28] obtain self-consistent solutions 
then? The reason is in the way of solving the equation of motion (129) in the nonlinear model. 
Given an appropriate initial form of F(r), one can evaluate {^'4>f^^ and {^^^TipY^^ by using 
the Kahana-Ripka plane- wave basis as long as the box size D and the maximum momentum 
kmax are finite. A new F{r) can then be obtained from 



f(r) = arctau| <^7^^>r |. (133) 

As krnax iucreases, both of < -^-^ >^ and < ipi^^Tip >^ tend to diverge. We numerically 
find that both quantifies increases at the same rate as k^ax increases so that the resultant 
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F(r) is quite insensitive to the value of k^ax for large enough kmax- This is the reason why 
stable sohton solutions could be found in the above mentioned single-subtraction Pauh-Villars 
regularization scheme. The existence of finite energy soliton could also be inferred from the 
derivative expansion analysis of the nonlinear lagrangian (123) with vanishing current quark 
masses. Nevertheless, one should keep in mind that it is not a completely satisfactory scheme 
in the sense that its predictions for some special quantities like the vacuum quark condensate 
contain divergences. For obtaining satisfactory answers also for these special quantities, the 
single-substr action Pauli-Villars scheme is not enough. We found that more sophisticated 
Pauli-Villars scheme with three substr action meets this requirement, and that its self-consistent 
solutions are only shghtly different from those of the naive single-subtraction scheme, except 
when discussing some special quantities as pointed out above. (This analysis will be reported 
elsewhere.) Considering the fact that the calculation of the structure functions are very time- 
consuming, we shall then use here the single-subtraction Pauli-Villars scheme, keeping in 
mind that some particular observables are out of the application of this regularization scheme. 
Finally, as for the nucleon mass, we prefer to using the theoretical soliton mass M^r ~ 1102 MeV 
rather than the physical mass, since it respects the energy-momentum sum rule at the energy 
scale of the model. 

For evaluating quark distribution functions at the O(O^), wc must perform infinite double 
sums over all the single-quark orbitals which are eigenstates of the static Hamiltonian H 
given by (31). As far as static nucleon observables, a numerical technique for carrying out 
such double sums was established in [5]. On the other hand, several new subtleties arising 
in the evaluation of quark distribution functions have been explained in [17]. In the actual 
numerical calculation, the expression of each physical quantity is divided into two pieces, i.e. 
the contribution of what we call the valence quark level (it is the lowest energy eigenstate of 
the static Dirac hamiltonian H, which emerges from the positive energy continuum) and that 
of the Dirac sea quarks (or the vacuum polarization contribution) as explained in [17]. The 
regularization is introduced into the latter part only. 

Now we start to show the results of our numerical calculation for polarized quark distribu- 
tion functions of the nucleon. Shown in Fig.l are the 0(0°) contributions to the isovector lon- 
gitudinally polarization distribution functions Au{x) — Ad{x) (solid curves) and Au{x) — Ad{x) 
(dashed curves), which was first calculated by Diakonov et al. [14]. Here Fig. 1(a) represents 
the contributions of the discrete valence quark level, while Fig.l (b) is the vacuum polarization 
contribution to the same quantities. Sum of these two contributions are shown in Fig. 1(c). 
As shown by Diakonov et al., the 0(il°) vacuum polarization contributions to Au{x) — Ad{x) 
and Au{x) — Ad{x) are fairly large. (The large and positive longitudinal polarization of the 
isovector combination of the antiquark distributions seems to be a characteristic prediction of 
the CQSM, which can in principle be tested by the improved phenomenological analyses of po- 
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larized parton distribution functions in the near future. More detailed discussion on this point 
will be given after finishing the evaluation of the O(Jl^) contribution to the same distribution 
function as well as that of the isoscalar longitudinally polarized distribution functions.) 

Next, we show in Fig. 2 the 0{Q^) contribution to the same distribution functions Au{x) — 
Ad{x) and Au{x) — Ad{x). Fig.2(a), 2(b) and 2(c) respectively stand for the 0{Q^) con- 
tributions of the discrete valence quark level, those of the Dirac sea quarks (or the vacuum 
polarization contributions), and their sums. One sees that the O(Jl^) contributions to the 
isovector longitudinally polarized distribution function are far from negligible as compared 
with the leading 0(fi°) contributions. This could be expected since the first moment of this 
distribution functions gives the isovector axial coupling constant of the nucleon. 



while we already know from the previous analyses that the 0{Q^) contribution to g)^ is 
large enough to resolve the longstanding qa problem in the hedgehog soliton model [29,30,24]. 
Adding this 0{Q^) contribution to the leading 0{Q^) contribution, wc obtain final answers 



answer for the isoscalar longitudinal distribution functions Au{x) + Ad{x) and Au{x) + Ad{x). 
Before showing those, we give in Fig. 3 the result for the 0{n^) contributions to the isoscalar 
longitudinally polarized distribution functions, Au{x) + Ad{x) and Au{x) + Ad{x). (We recall 
that there is no 0{DP) contribution to these distribution functions.) Fig.3(a), 3(b) and 3(c) 
respectively stand for the contributions of the discrete valence quark level, those of the vacuum 
polarization contributions, and their sums. One sees that the vacuum polarization contribu- 
tions to the distribution functions Au{x) + Ad{x) and Au{x) + Ad{x) are much smaller than 
those of the corresponding isovector distributions Au{x) — Ad{x) and Au{x) — Ad{x). Now we 
show in Fig.4 the final answers for Au{x) — Ad{x) and Au{x) —Ad{x), which are the sums of the 
0{DP) and O(Q^) contributions, in comparison with the final answers for Au{x) + Ad{x) and 
Au{x) + Ad{x) arising from the 0{Q^) terms alone. We observe quite a big difference between 
the isovector distributions and the isoscalar one. The overall magnitude of Au{x) + Ad{x) is 
much smaller than that of Au{x) — Ad{x), which denotes that li-quark is positively polarized, 
while the d-quark is negatively polarized to the direction of proton spin. 

At this stage, it may be interesting to compare our theoretical predictions for the lon- 
gitudinally polarized quark distribution functions with some of the semi-phenomenological 
parametrization. The parametrization given by Gliick, Reya, Stratmann and Vogelsang is es- 
pecially convenient for the purpose of handy comparison [31], since the normalization point 
{Qlnit — 0.34 GeV^) of their parametrization is fairly close to the energy scale of our effec- 
tive quark model {Mpy ~ 0.32 GcV^. Fig.5 shows this comparison. The filled squares in 
Fig. 5(a) and Fig.5(b) stand for the GRSV parametrizations for the quark distribution func- 
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for Au{x) — Ad{x) and Au{x) — Ad{x), which will be shown later together with the final 
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tions x{Au{x) + Au{x) + Ad{x) + Ad{x)) and x{Au{x) + Au{x) —Ad{x) —Ad{x)), respectively. 
Of the two theoretical curves in each figure, the sohd curve is the answer of the present cal- 
culation, whereas the dashed curve is obtained by using the old treatment used in [17], which 
amounts to dropping some of the nonlocality effects in time. One observes that the nonlo- 
cality corrections newly introduced in the present analysis are quite important especially for 
the isoscalar distribution x{Au{x) + Au{x) + Ad{x) + Ad{x)), while it is less important for 
the isovector distribution x{Au{x) + Au{x) — Ad{x) — Ad{x)). (This is probably because the 
nonlocahty correctios appearing at the 0{Q^) are masked by the dominant contribution 
in the case of isovector polarized distribution functions.) By comparing the two theoretical 
curves for x{Au{x) + Au{x) + Ad{x) + Ad{x)) with the corresponding GRSV parametrization, 
one finds that the new treatment leads to a better agreement. Especially impressive is that the 
new treatment reproduces the negative sign of the GRSV distribution function in the smaller 
x region, although one should not forget the fact that the GRSV parametrizations are not 
experimental data themselves. We point out that the most important factor leading to this 
qualitative difference between the old and new treatments of the quark distribution functions 
is the nonlocality correction arising from the second term of (38), i.e. the proper account of 
nonlocality in time of the operator A^{0)OaA{zo). 

Turning back to Fig. 4, let us inspect the theoretical predictions for the antiquark distribu- 
tions in more detail. An interesting feature is that, in most region of x, Au{x) — Ad{x) > 
and Au{x) + Ad{x) < with the relation \ Au{x) — Ad{x)\ ^ \ Au{x) + Ad{x)\. This denotes 
that d is strongly polarized in the opposite direction to the proton spin, while u is weakly 
polarized in the same direction to the proton spin. This appears to be a prominent prediction 
of the CQSM, which is worthy of special mention. In fact, it sharply contradicts the assump- 
tion of SU(2) symmetric sea quark polarization Au{x) = Ad{x), which is frequently used in 
semi-phenomenological analyses of parton distributions. The isospin symmetric polarization 
is also assumed in the analysis by Gliick, Reya, Stirling and Vogt [31]. We compare in Fig. 6 
our prediction for the xAu{x) and xAd{x) with the GRSV parametrization, which assumes 
that xAu{x) = xAd{x) (= xAq{x)). Naturally, one finds qualitative difference between the 
theoretical distributions and the GRSV parametrization. Still, it is interesting to see that 
the average of the two theoretical distributions xAu{x) and xAd{x) is not extremely different 
from the corresponding GRSV parametrization xAq{x). As for the unpolarized distribution 
functions, the breakdown of the assumption of SU(2) symmetric sea has already been con- 
firmed by the NMC measurement [16]. By the same token, there is no compelling reason to 
believe that the spin dependent antiquark (sea quark) distributions are isospin symmetric. In 
fact, our previous analyses based on the same model shows that the isospin asymmetry of the 
unpolarized sea quark distributions can be explained very naturally as combined effects of two 
ingredients, i.e. the apparently existing flavor asymmetry of valence quark numbers in the 
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nucleon and the spontaneous chiral symmetry breaking of QCD vacuum [17-19]. It is just the 



The above-mentioned fairly big difference between the isovector and isoscalar longitudinally 
polarized distribution functions manifests itself also in their first moments, i.e. the isovector 
and isoscalar axial charges given as 



The resultant large isovector axial charge and small isoscalar (flaver-singlet) one seem to be 
qualitatively consistent with the observation. Especially interesting here is the flavor- singlet 
axial charge identified with the quark spin content of the nucleon. In the context of the 
CQSM, this quantity was first investigated in [5] with use of the self-consistent soliton solution 
obtained in the proper-time regularization scheme. The value of g^^^ = (S3) obtained there 
ranges from 0.4 ~ 0.5 corresponding to the variation of the dynamical quark mass M from 
425 MeV to 375 MeV. One may notice that the value g^^^ ~ 0.35 obtained in the present 
calculation is a little smaller than the previous one. The cause of this difference can be 
traced back to the qualitative change of the self-consistent soliton solution obtained in the 
new regularization scheme. As a general trend, the Pauli-Villars regularization scheme cuts off 
high momentum components more weakly than the energy-cutoff scheme like the proper-time 
one, thereby leading to soliton solutions with stronger distortion. Incidentally, owing to the 
nucleon spin sum rule (L3) + ^(Ss) = ^ proved in [5], the rest of the nucleon spin is carried 
by the orbital angular momentum of the effective quark fields. (Naturally, this is true only at 
low corresponding to the energy scale of our effective model. It will be shown later that 
an increasing portion of the nucleon spin is carried by gluons as increases.) A soliton with 
stronger distortion gives larger orbital angular momentum, and consequently smaller quark 
spin fraction [5]. 

The characteristic feature of the above theoretical prediction, i.e. larger isovector charge 
and smaller isoscalar one seems also consistent with the idea of Nc counting or 1/A^c expansion 
of QCD. For understanding it, we just recall the fact that the collective angular velocity 
Q scales as 1/A^c, so that the leading contributions to the isovector and isoscalar polarized 
distribution functions are respectively of the 0{N^) and 0{N^). The detailed comparison of 
the theoretical first moments with the corresponding experimental data will be given later 
after taking account of the scale dependence of them. 

Now we show the results of our numerical calculation for transversity distributions. Fig. 7 
shows the 0(0°) contributions to the isovector transversity distribution functions 6u{x) — Sd[x) 
(solid curves) and Su{x) — Sd{x) (dashed curves). Here Fig.7(a) stands for the contributions 



same mechanism that is responsible for the opposite longitudinal polarization of the u and d 
quarks. 
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of the discrete valence level, while 7(b) is the vacuum polarization contributions to the same 
quantities. The sum of these two contributions are shown in Fig. 7(c). One finds that the vac- 
uum polarization contributions to these distribution functions are insignificant. Next, we show 
in Fig. 8 the O(Jl^) contributions to the same transversity distribution functions 5u{x) — 5d{x) 
and 5u{x) — 5d{x). Also for these O(Q^) terms, the vacuum polarization contributions are 
very small as compared with the contributions of the discrete valence quark level. However, 
we emphasized that the valence level contribution at the O(Q^) is far from small as compared 
with the leading 0(fi°) contributions, and should not be discarded. Shown in Fig. 9 are the 
theoretical isoscalar transversity distributions resulting at the O(Q^). One sees that vacuum 
polarization contributions to the isoscalar transversity distributions is also rather small. The 
final predictions of the CQSM for 5u{x) — 5d{x) and 5u{x) — Sd{x), which are the sums of the 
0{Q^) and 0{Q^) contributions are shown in Fig. 10, in comparison with the final answer for 
Su{x) + 6d{x) and 6u{x) + 6d{x) arising from 0{fl^) terms. One again sees that the magnitudes 
of the isoscalar distributions are much smaller than those of the isovcctor distributions. Re- 
member the similar observation made for the longitudinally polarized distribution functions. 
To see it in more detail, we find that the ratio of the isoscalar to isovector distribution are 
much smaller for the longitudinally polarized distribution than for the transversity one. We 
shall come back later to this point when discussing the corresponding first moments of these 
spin dependent quark distribution functions. 

Roughly speaking, the quark distribution functions evaluated here corresponds to the en- 
ergy scale of the order of the Pauli-Villars cutoff mass Mpy ^ 0.56 GcV. The evolution 
must be taken into account in some way before comparing them with the observed nucleon 
structure functions at high Q^. Recently, Saga group provided a Fortran program, which gives 
numerical solution of DGLAP evolution equations at the next-to-leading order (NLO) for the 
polarized as well as unpolarized structure functions of the nucleon [32-34] . We shall make use 
of their Fortran programs to evaluate the polarized distribution functions at large [33,34]. 
The question here is what value we should take for the initial energy scale of this Q'^ evolution. 
Since the use of perturbative QCD below 1 GeV is anyhow questionable, one may take this 
initial energy scale Q'^^u ^ adjustable parameter, which would be fixed by adjusting the 
observed structure functions at high energy region. Here we have tried to see the effect of 
variation of Qi^u ^ small range of around the model energy scale of Mpy ~ (0.56 GeV)^. 
The value Qinit — (0-5 GeV)^ = 0.25 GeV^ obtained from this analysis will be used throughout 
the following investigation. Before showing the results of evolution, we want to make a 
short comment. One notices from the figures given so far, the distribution functions evaluated 
in our effective model have unphysical tails beyond x > 1, although they are not so significant. 
These unphysical tails of the theoretical distribution functions come from an approximate 
nature of our treatment of the soliton center-of-motion (as well as the collective rotational 
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motion), which is essentially nonrelativistic. A simple procedure to remedy this defect was 
proposed by Jaffe based on the (1 + 1) dimensional bag model [35] and recently reinvesti- 
gated by Gamberg et al. within the context of the N JL soliton model [36] . According to the 
latter authors [36] , the effect of Lorentz contraction can effectively be taken into account by 
first evaluating the distribution functions in the soliton rest frame (as we are doing here) and 
then by using a simple analytical transformation that preserves first moments of distribution 
functions, as far as the 0(Q°) contributions to the distribution functions are concerned. Such 
a simple relation may not be expected however if we consider the rotational motion of the 
soliton, which are anyhow three dimensional. In fact, a comparison with the corresponding 
phenomenological distribution functions seems to indicate that the above procedure based on 
the (1 + 1) dimensional dynamics tends to overestimate the effect of Lorentz contraction. In 
the present investigation, we therefore decided not to use their procedure. Still we want distri- 
bution functions which vanish outside the range < a; < 1 so that we can use the Q^-evolution 
Fortran program provided by Saga group [33,34]. Since the unphysical tails of our theoretical 
distributions are rather small in magnitude, we are to use a simple cutoff procedure as follows. 
That is, we obtain modified distribution functions, which can be used as input distributions 
of the above Fortran program, from the original theoretical distribution functions by multi- 
plying the x-dcpcndcnt cutoff factor (1 — x^^). (This special cutoff factor is invented from the 
requirement that only the tails of the distribution functions are modified.) Fig. 11 illustrates 
the effect of this tentative cutoff procedure. The solid curve here is the theoretical distribution 
function Au{x) — Ad{x) given as a sum of the 0(f2°) and 0{Q^) contributions. We point out 
that this distribution function is the worst case in the sense that the tail beyond a; = 1 is most 
significant as compared with the other distributions. The dashed curve in the same figure is 
obtained by using the above cutoff procedure. One sees that it leaves the distribution function 
for X < 0.7 almost intact. Naturally, this cutoff procedure alters the values of integrals of the 
distribution functions, i.e. the first moments. However, it turns out that the reduction is less 
than 2 % even in the above worst case. We therefore expect that the tentative nature of the 
above procedure hardly affects the following qualitative analyses of scale dependence of the 
quark distribution functions. 

For the sake of comparison, we have carried out a similar evolution procedure also for the 
initial distributions given by the MIT bag model. The distribution functions of the (naive) 
MIT bag model are already known and they are given analytically as follows [22] . The isoscalar 
longitudinally polarized distribution functions is given by 




(137) 
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whereas the isoscalar transversity distribution functions is given as 
8u[x)+M[x) = — — 7T-2Z-T / dyy to{uJi,y) 



+ 2io(a;i,t/)tiKy)(^) + t?(a;i,y)(^r]}. (138) 



y y 

On the other hand, the isovector distribution functions are simply related to the isoscalar ones 
as 

Au{x) - Ad{x) = ^ [Au{x) + Ad{x) ], (139) 

6u{x)-dd{x) = ^ [du{x) + Sd{x)]. (140) 

In eqs. (137) and (138), Mjv and R respectively stard for the nucleon mass and the bag radius, 
while cun is the nth root of the bag eigenvalue equation as 

tana;„ = (141) 

UJn — i- 

and ymin = \xMj^R — uJi\. The function t/(a;„, y) is defined by 

ti{ujn,y)= I ji{uujn)ji{uy)u'^du. (142) 
Jo 

The bag radius R is only one free parameter of this simple model. In the numerical calculation, 
we adopt the value used by Jaffe and Ji [22], i.e. 

MnR = 4:.Qu;i, (143) 

where loi ~ 2.043 is the lowest (dimensionless) eigenvalue of the bag equation. 

To get a rough idea about the scale dependence, we show in Fig.l2 and Fig.l3 the theoretical 
polarized quark distribution functions before and after Q^-evolution. Here Au{x) and 5u{x) 
in Fig. 12 (a) respectively stand for the longitudinal and transversity distributions for li-quark. 
In our model, the difference between the two distributions are sizable even at the initial low 
energy scale. A comparison with the existing and yet-to-be-obtained high energy data must 
be done with care, since the way of evolution of these two distributions are pretty different 
and the deference between the two becomes larger and larger as Q"^ increases. A general 
trend is a rapid growth of small x component of the longitudinally polarized distribution due 
to the coupling with gluons. A similar tendency is also observed for the corresponding d- 
quark distributions shown in Fig. 12(b). We can also give some predictions for the polarized 
antiquark distribution functions. As one can sec in Fig. 13, even the signs are different for 
the longitudinal and transversity distributions. (This is the case for both of u and d quarks.) 
The twist-2 spin dependent distribution functions were calculated by several authors based 
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on various effective models of baryons [38,39,12]. As for the polarized quark distribution 
functions, the predictions of various models give more or less similar shape of distributions 
assuming that they take account of the dominant nature of the valence quark contribution as 
well as the effects of pion cloud in some effective way. The situation is quite different for the 
polarized antiquark distributions. The transversity distribution functions for the antiquarks 
have, for instance, been evaluated by Barone et al. within the chiral chromodielectric model 
[38]. Comparing their predictions for 5u{x) and 5d{x) with ours shown in Fig. 13, we find that 
their model gives 6u{x) > 0, while ours does Su{x) < 0. The shapes of Su{x) and Sd{x) are 
also quite different in both models. In consideration of the fact that the polarized antiquark 
distributions are quite sensitive to the detailed dynamics of the model, it is very important to 
get precise phenomenological information for them. 

Next we show in Fig. 14(a) the theoretical predictions for the proton structure function 
gi{x, Q^) at = 5 GeV^ in comparison with the corresponding experimental data given by 
E143 collaboration [40]. The theoretical curves are obtained as follows. Starting with the initial 
distributions Au{x) + Ad{x) , Au{x) + Ad{x) and Au{x) — Ad{x) , Au{x) — Ad{x) or equivalcntly 
Au{x),Au{x) and Ad{x),Ad{x) given at QLit — 0.25 GeV^ (we assume As{x) = As{x) = 
and Ag{x) = at this energy scale), we solve the NLO evolution equation to obtain the dis- 
tribution functions at = 5 GcV^. These distribution functions are then convoluted with 
the relevant quark and gluon coefficient functions at the NLO within the framework of pcrtur- 
bative QCD. These procedures have been carried out for the initial distribution given by the 
CQSM and also by the MIT bag model. The solid and dashed curves in Fig. 14(a) respectively 
stand for the prediction of the CQSM and that of the MIT bag model. A remarkable feature of 
the CQSM as compared with the MIT bag model is the enhancement of the structure function 
at small x region, i.e. large sea quark components. One also observes that a clear peak of 
gi{x, Q^) around x ~ 0.3 predicted by the MIT bag model (a relativistic valence quark model) 
is not seen in the experimental structure function. On the other hand, one can say that the 
prediction of the CQSM reproduces qualitative feature of the observed structure function in 
the whole range of x. Fig. 14(b) shows the theoretical prediction of the CQSM (sohd curve) and 
that of the MIT bag model (dashed curve) for the neutron spin structure function gi{x, Q^) in 
comparison with the E154 data [41]. One clearly sees that the neutron spin structure function 
gi{x, Q^) predicted by the MIT bag model is neghgibly small in magnitude even after evolution. 
We recall that at the initial energy scale the naive MIT bag model predict gi{x) — 0, which 
is a necessary consequence of a model that does not properly incorporate chiral symmetry. 
On the other hand, the prediction of the CQSM for gi{x, Q^) is seen to be large and negative 
especially in the small x region in good agreement with the experimental observation. Then, 
this agreement may be regarded as a manifestation of the importance of chiral symmetry in 
the physics of high-energy deep-inelastic scattering. 
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As is widely known, the simplest but the most important quantities characterizing the 
quark distribution functions are the associated first moments. Here we are interested in the 
first moments of the longitudinally polarized distribution functions and of the transversity 
ones, which are respectively called the axial and tensor charges defined as 

gf = {[^u(x) - Ad(x)] + [Au(x) - Ad(x)]} dx, (144) 
Jo 

gf = / {[Au(x) + Ad(x)] + [Au(x) + Ad(x)]} dx, (145) 

^0 

= / {[5u(x)-5d(x)]-[5u(x)-5d(x)]}dx, (146) 
Jo 

gP = / {[5uix)+5dix)\-[5uix)+5dix)]}dx. (147) 

Before discussing the prediction of the CQSM for these quantities, it may be instructive to 
remember some basic properties of those. (We recall that the first calculation of the tensor 
charge in the CQSM was given in [42] .) As emphasized by Jaffe and Ji [22] , there is a remark- 
able difference between the axial and tensor charges originating from the charge conjugation 
properties of the relevant operators. For each flavor, the tensor charge counts the number 
of valence quarks (quarks minus antiquarks) of opposite transversity. Consequently, the sea 
quarks do not contribute to the tensor charge. (This does not necessarily means vanishing 
transverse polarization of antiquarks, however.) On the other hand, the axial charge counts 
the number of quarks plus antiquarks of opposite helicity. In fact, by rewriting (146) and (147) 
as 

gf^ = {[Au{x) - Ad{x)] - [Au{x) - Ad{x)]] dx + 2 f\Au{x) - Ad{x)] dx, (148) 
^0 Jo 

= / {[Au(x) + Ad(x)] - [Au(x) + Ad(x)]} dx + 2 f [Au(x) + Ad(x)] dx, (149) 
JO Jo 

the first and the second terms of the above equation can respectively be interpreted as valence 
and sea quark contributions in the parton model. Since the sea quark degrees of freedom is 
absent in the nonrelativistic framework, the difference between the axial and tensor charges is 
purely relativistic. Still, one must clearly distinguish two types of relativistic effect. The one 
is dynamical effects, which generate sea quark polarization. The other is kinematical effects, 
which make a difference between the axial and tensor charges even though the sea quark degrees 
of freedom arc totally neglected. The existence of this latter effect can readily be convinced by 
comparing the prediction of two "valence quark models", i.e. the nonrelativistic (constituent) 
quark model and the MIT bag model. In fact, the non-relativistic quark model predicts 

9f - # - I (150) 

- gP = 1, (151) 
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while the prediction of the MIT bag model is given by 

g^T = I- j{f~\9')r'dr, gP^l- j{f + \9')r'dr, (152) 

= i-Jif-lg')r'dr, gP = i-J{f + lg')r'dr, (153) 

where / and g are upper and lower components of the lowest energy quark wave functions. 
For a typical bag radius R ~ 4.0a;i/Mjv, which was used before, this gives 



gf ~ 1.06, g'^' ~ 1.34, 

0.84. 



,(3) 



gA 



0.64, # 



(154) 
(155) 



As is obvious from eqs. (152) and (153), the splittings of the axial and tensor charges are due 
to the different sign of the lower component (p-wave) contributions [22]. One should however 
notice that there is one interesting feature shared by both of the nonrelativistic quark model 
and the MIT bag model. The predictions of the both models for the ratio of the isoscalar to 
isovector axial charges as well as the ratio of the isoscalar to isovector tensor charges are just 
the same : 

gf/gf = gP/g? = 3/5. (i56) 

Although there is no experimental information yet for the tensor charges, the above prediction 
for the ratio of the two axial charges obviously contradicts the EMC observation. Now we 



Table 1: The theoretical predictions for the isovector and isoscalar axial charges as well as the 
corresponding tesnsor charges. The predictions of the MIT bag model and those of the lattice 
QCD [43] are also shown together with some experimental data [44,45]. 





CQSM 


MIT-bag 


Lattice QCD [43] 


Experiment 




1.41 


1.06 


0.99 


1.254 ± 0.006 [44] 
(Q^-indep.) 




0.35 


0.64 


0.18 


0.31 ± 0.07 [45] 
(g2 = 10 GeV") 


(3) 

gV 


1.22 


1.34 


1.07 




g'? 


0.56 


0.80 


0.56 




g'^'lgT 


0.25 


0.60 


0.18 


0.24 


g^^^/g? 


0.46 


0.60 


0.52 
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shall argue that the above prediction may be interpreted as showing the limitation of simple 
valence quark models, which fail to properly incorporate chiral symmetry of QCD. To convince 
it, the predictions of the NRQM and the MIT bag model are compared with those of the 
CQSM in table 1, which maximally incorporate chiral symmetry. For the sake of reference, the 
predictions of the lattice QCD are also shown [43]. (Here we have omitted the errors of the 
lattice QCD calculation, for simplicity.) We first point out that the predictions of the CQSM 
for the above ratios, i.e. 

gf/gf 0.25, 0.46, (157) 

strongly deviate from the above predictions of the two valence quark models. What is re- 
markable here is that the CQSM predicts very small isoscalar axial charge in consistent with 
the EMC observation. (More meaningful comparison should be made after taking account 
of the scale dependence of this quantity.) Its prediction for the isovector axial charge is also 
qualitatively consistent with the experimental value determined from the neutron beta decay. 
(The deviation from the experimental value is only about 11 %.) The lattice gauge theory 
also predicts a very small isoscalar axial charge g^^^ ~ 0.18. However, this prediction may not 
be taken as a final one since it largely underestimates the isovector axial charge. At any rate, 
one can observe qualitative similarities between the predictions of the CQSM and those of the 
lattice QCD. Both predicts quite a small number for the ratio of the isoscalar to isovector 
axial charges as compared with the prediction g^^ / Qa' — 0-6 of the NRQM or the MIT bag 
model. On the other hand, the predictions of both models for the ratio of the isoscalar to 
isovector charges is not extremely different from the prediction g^^ / g^^ = 0.6 of the latter 
valence quark models. In our opinion, the observed deviation from the valence quark picture 
indicates an importance of chiral symmetry as a generator of "dynamical sea quark effect" , and 
the predicted feature is expected to be confirmed by future measurements of tensor charges. 

To compare the theoretical first moments of the spin distribution functions with the existing 
data for the longitudinal case and with yet-to-be-observed ones for the transversity case, we 
must take account of the scale dependence of the relevant moments. As is well-known, the 
first moment of the isovector longitudinal distribution functions, i.e. the isovector axial charge 
is scale independent, i.e. it does not evolve : g^A^Q"^) — 9A\Qinit)- This is due to the 
conservation of the fiavor nonsinglet axial- vector current [51]. This is not generally the case 
for the fiavor singlet (isoscalar) axial charge owing to the so-called axial anomaly of QCD 
[46,47]. (Still, one can take a scheme called the chiral invariant factorization scheme in which 
the flavor singlet axial charge is independent of [48]. Here, we take more standard gauge 
invariant factorization scheme [49] .) In the singlet sector, the nth moments of the longitudinally 
polarized distribution functions are coupled with the corresponding gluon contributions. The 
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evolution of these nth moments is governed by the anomalous dimension matrix 



igq 199 



where 7^°^ and 7^^^ are 1- and 2-loop contributions to the anomalous dimensions. An analytic 
solution to this coupled evolution equation of the NLO is given in the matrix form [50,51] : 



r"(^3') = I 71:1 1 , (159) 



+ 



pn _ 1 Cy.s{Qinit) CX.sjQ ) pn^npr. 



' asm ' 

^ Olg (Qinit) , 



2/3o 



Air 



pn^n pn 



47r 



47r WiQ^ 



init) 



(160) 



Here as{Q^) is the QCD running coupling constant at the next-to- leading order with MS 
scheme, /So and j3i are the 1- and 2-loop QCD beta functions, respectively, and 



^ ^(l)n _ !^^(0)n_ 

P" are 2x2 projection matrices defined by 
with i being a 2 x 2 unit matrix and with 



'99 



(0)n {0)n 



(161) 



(162) 



(163) 



the eigenvalues of the 1— loop anomalous dimension matrix 7^'^)". Since the necessary anoma- 
lous dimension matrices are all given in [51], it is easy to calculate the evolution of the first 
moment of the flavor singlet longitudinally polarized distribution functions, i.e. the isosinglet 
axial-charge. 

Because of its chiral-odd nature, the moments of the transversity distributions do not 
couple with gluons, irrespective of the flavor quantum numbers, which especially means that 
isovector and isoscalar tensor charges follow the same evolution equation. The anomalous 
dimension of the transversity distribution at the leading 1—loop order was first given by Artru 
and Mekhfi [52], while the corresponding 2— loop contributions have recently been given by 
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three groups independently [53-55] . Once the relevant anomalous dimensions are known, it is 
easy to obtain an analytical solution of the NLO evolution equation for the nth moment of 
transversity distribution. Here, we use the form given by Hayashigaki et al. [54] as 

(l)n (0)n 

where the relevant anomalous dimensions 7;°''" and 7}^^''" arc all given in [51]. Fig. 15 show the 
calculated dependence of the axial and tensor charges. For obtaining it, we start with the 
theoretical first moments given at the initial energy scale Q^^u ~ 0.25GeV : 

gfiQlnt) = 1-41 (165) 

gfiQlu) ^ AE(QLt) = 0.35 (166) 

AG{Ql,,) = (167) 

9?\Qlu) = 1-22 (168) 

gPiQLt) = 0.56 (169) 

One sees that the dependence of the flavor singlet axial charge is very small (it is almost 
constant except in the very low region). A characteristic prediction of the CQSM for the 
axial charges, i.e. large isovector charge and small isoscalar charge appears to be qualitatively 
consistent with the corresponding experimental data at the relevant energy scale. As was 
pointed out by many authors [38,52-55], the dependence of the tensor charges are sizably 
large. Although there is no experimental information for these latter quantities, this depen- 
dence must be taken seriously when comparing the theoretical prediction of low energy models 
with future experimental date. (Note however that that the ratio g^^ / g^^ is independent.) 

Because of the coupling between the flavor singlet axial charge (the longitudinal quark 
polarization) and the gluon polarization in the evolution equation, non-zero gluon polarization 
appears at high even if we have assumed AG = at the initial energy scale of Ql — 
0.25 GeV^. We show in Fig. 16 the Q'^ evolution of AG in comparison with that of AE = 
g^^\ One sees that the gluon polarization rapidly grows with increasing Q^. Already at 

~ 2GeV^, AG is seen to be larger than AE. As explained in [49], the growth of the 
gluon polarization with can be traced back to the positive sign of the anomalous dimension 
7^°)-'^ at the leading order (7^^^^ = 2). The positivity of this quantity means that a polarized 
quark is preferred to radiate a gluon with helicity parallel to the quark polarization. Since 
the net quark spin component in the proton is positive, it follows that AG > at least for 
the gluons perturbatively emitted from quarks [49]. It is hoped that the direct information 
on Ag[x, Q'^) from the di-jet asymmetry analyses at HERA in conjunction with the precise 
NLO analyses of gi{x, Q^) will soon provide us with an accurate determination of the polarized 
gluon distribution as well as its flrst moment [56] . 
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5 Summary 



In summary, we have shown that the CQSM naturally explains qualitative behavior of the 
experimentally measured longitudinally polarized structure functions of the proton and the 
neutron. As was shown in our previous papers, the model also reproduces an excess of d sea 
over the u sea in the proton very naturally [17-19]. Furthermore, it predicts quahtative differ- 
ence between the transversity distribution functions and longitudinally polarized distribution 
functions. For example, in simple valence quark models like the NRQM or the MIT bag model, 
the ratios of the isoscalar to isovector charges are just the same for both of the axial charges 
and the tensor charges. On the contrary, in the CQSM or in the lattice gauge theory, this 
ratio turns out to be much smaller for the axial charges than for the tensor charges. In our 
viewpoint, what makes this difference is "dynamical sea quark effects" dictated by the spon- 
taneous chiral symmetry breaking of the QCD vacuum. Another noteworthy prediction of the 
CQSM is the opposite (spin) polarization of the u and d sea quarks, thereby indicating SU(2) 
asymmetric sea quark polarization. These observations then indicates that nonperturbative 
QCD dynamics due to the spontaneous chiral symmetry breaking would survive and manifest 
itself in the isospin (or flavor) dependence of high energy spin observables, especially in that 
of the polarized (as well as unpolarized) antiquark distribution functions. 
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Figure caption 



Fig. 1. The 0{0P) contributions to the isovector longitudinally polarized dis- 
tribution functions Au{x) — Ad{x) (solid curves) and Au{x) — Ad{x) 
(dashed curves). Here the three figures (a), (b), and (c) correspond to 
the contributions of the discrete valence level, those of the Dirac sea 
quarks, and their sums, respectively. 



Fig. 2. The 0{fl^) contributions to the isovector longitudinally polarized dis- 
tribution functions Au{x) — Ad{x) (solid curves) and Au{x) — Ad{x) 
(dashed curves). The meaning of the three figures (a), (b), and (c) is the 
same as in Fig.l. 



Fig. 3. The O(Q^) contributions to the isoscalar longitudinally polarized distri- 
bution functions Au{x)+Ad{x) (solid curves) and A'u(a;)-|-A(i(a;) (dashed 
curves). The meaning of the three figures (a), (b), and (c) is the same as 
in Fig.l. 



Fig. 4. The final predictions of the CQSM for the longitudinally polarized dis- 
tribution functions Au{x) — Ad{x) and Au{x) — A(i(x') given as the sums 
of the 0{Vf) and 0{p}) contributions, in comparison with those for the 
isoscalar longitudinally polarized distribution functions Au{x) -\- Ad{x) 
and Au{x) + Ad{x) comong from the 0{Q}) terms. 



Fig. 5. The theoretical predictions for the longitudinally polarized distribu- 
tion functions, x{Au{x) + Au{x) -\- Ad{x) + Ad{x)) and x{Au{x) -\- 
Au{x) — Ad{x) — Ad{x)), are compared with the corresponding semi- 
phenomenological parametrization given by Gliick, Reya, Stratmann and 
Vogelsang [31]. Of the two theoretical curves in each figure, the solid 
curve is the answer of the present calculation, whereas the dashed curve 
is obtained by using the old treatment used in [17], which amounts to 
dropping some of the nonlocality effects in time. 
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Fig. 6. The predictions of the CQSM for the polarized antiquark distribu- 
tions xAu{x) and xAd{x) are compared with the corresponding GRSV 
parametrization, which assumes SU(2) symmetric sea quark polarization, 
xAu{x) — xAd{x) (= xAq{x)). 

Fig. 7. The 0(f2°) contributions to the isovector transversity distribution func- 
tions Su{x) — 5d{x) (solid curves) and 5u{x) — Sd{x) (dashed curves). 
Here the three figures (a), (b), and (c) correspond to the contributions of 
the discrete valence level, those of the Dirac sea quarks, and their sums, 
respectively. 

Fig. 8. The 0{Q^) contributions to the isovector transversity distribution func- 
tions 6u{x)—5d{x) (sohd curves) and 5u{x) — Sd{x) (dashed curves). The 
meaning of the three figures (a), (b), and (c) is the same as in Fig.l. 

Fig. 9. The 0(f2^) contributions to the isoscalar transversity distribution iunc- 
tions 5u{x) + Sd{x) (solid curves) a,nd Su{x) + 5 d{x) (dashed curves) . The 
meaning of the three figures (a), (b), and (c) is the same as in Fig.l. 

Fig. 10. The final predictions of the CQSM for the transversity distribution func- 
tions 5u{x) —Sd{x) and Su{x) — Sd{x) given as the sums of the 0(0") and 
0{fl^) contributions, in comparison with those for the isoscalar transver- 
sity distribution functions 5u{x) + 5d{x) and 5u{x) + 5d{x) coming from 
the 0{n^) terms. 

Fig. 11. The solid curve represents the theoretical distribution functions Au{x) — 
Ad{x), whereas the dashed curve is a modified one obtained from it by 
multiplying a x-dependent cutoff factor (1 — x^°). 

Fig. 12. The theoretical predictions for the twist-2 spin dependent quark distri- 
bution functions before and after Q^-evolution. Here Au{x) and 6u{x) 
(in (a)) respectively stand for the longitudinal and transversity distribu- 
tions of li-quark, while Ad{x) and 5d{x) (in (b)) are the corresponding 
quantities for d-quark. 
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Fig. 13. The theoretical predictions for the twist- 2 spin dependent antiquark dis- 
tribution functions before and after Q^-evolution. Here Au{x) and 5u{x) 
(in (a)) respectively stand for the longitudinal and transversity distribu- 
tions of ■u-quark, while Ad{x) and 5d{x) (in (b)) are the corresponding 
quantities for d-quark. 

Fig. 14. The theoretical predictions for the proton and neutron spin structure 
functions g^{x,Q'^) and (a;, Q^) at = 4GeV^ in comparison with 
the corresponding SLAG data. The sohd and dashed curves in (a) re- 
spectively stand for the prediction of the CQSM and that of the naive 
MIT bag model for (x, Q'^) , whereas the black circles are the E143 data. 
The corresponding theoretical predictions for the gi{x, Q^) are shown in 
(b) together with the E154 data. 

Fig. 15. The scale dependence of the axial and tensor charges. The evolution 
equations at the next-to-leading order are solved under the initial con- 
ditions gf{Qj^^) = 1.41, gf\Qj^^) ^ AE(QLJ = 0.35, #(gLt) = 
1.22, #(QL,) = 0.56, and AG'(gLJ = at g^t = 0.25 GeV^ 

Fig. 16. The scale dependence of the flavor singlet axial charge (or the quark 
polarization) and the gluon polarization. The initial conditions for the 
evolution equation is the same as given in Fig. 15. 
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